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The nonlinear partial differential equations describing the spin dynamics of Heisen- 
berg ferro and antiferromagnet are studied by Lie transformation group method. The 
generators of the admitted variational Lie symmetry groups are derived and conser- 
vation laws for the conserved currents are found via Noether's theorem. 
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I. INTRODUCTION. 

As it has been demonstrated previously in a number of works the symmetry methods are 
very efficient tool in application to differential equations in physics. A subject of an especial 
interest is a study of invariance properties (symmetries) of the equations with respect to 
local Lie groups point transformations of dependent and independent variables. The de- 
tailed discussion and applications of the symmetries of differential equations and variational 
problems may be found in Ref.^ and 3 -. 

The non-linear dynamics of n-dimensional (anti)ferromagnet is described by non-linear 
second-order partial differential equations (PDE). In our previous work it has been shown 
that the system of differential equations are the Euler-Lagrange equations of a certain action 
functional and the corresponding long-wave action has been constructed via Vainberg's the- 
orem (see also Appendix). It is well known that a part of one parameter symmetry groups 
of these equations turns out to be their variational symmetries as well, i.e. a symmetries of 
the associated action functional. According to Noether theorem^ such an invariance of the 
elementary action is a necessary and sufficient condition of an existence of conservation laws 
for the smooth solutions of the initial PDE. 

The layout of the paper is as follows. By using the Noether identity we define the 
variational symmetries among the point Lie symmetries of the initial differential equations 
of a 2D-dimensional ferromagnet. After they are deduced we establish explicit formulae for 
the conserved densities and density currents involved in the conservation laws both for a 
Heisenberg ferromagnet and antiferro magnet. Apparently, the analysis presented may be 
easily generalized for an arbitrary dimension n > 3. 

II. DYNAMICAL INVARIANTS OF 2D FERROMAGNET. 

The non-linear partial differential equations in 3 independent variables (x 1 ,^ 2 ,^ 3 ) = (£, 
r) and two dependent variables u l = 9(t, r) and u 2 = Lp(t, r) describing the dynamics of 
Heisenberg ferromagnet in continuum approximation can be written as follows 4 

F 1 = —S sin 6if t +/3(a9- cos 9 sin 9 (Vtp) ^ = 0, 

F 2 = S sin 99 1 + /3 (2 cos 9 sin 9 (vWp) + sin 2 9A<p\ =0, (1) 
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where (3 = ^4-. Here and throughout the following notations are taken: a have the range 1, 
2 meaning 9 and (p, correspondingly, the usual summation convention over a repeated index 
is employed, 9j u ..j h (or <fij l ...j k ) denote the k-th. order partial derivatives of the dependent 
variables 

e =e dH 

(J) 31 - Jk dxndxj*...dx**' 
The J denotes a multi-index J = (ji,j2, ■ ■ ■ ,jk) with jf. = 1, . . . , 3 (k > 0) pointing to the 
independent variables with respect to which one differentiate. 

Consider a local one-parameter Lie group of point transformations acting on a space 
Q of the independent and dependent variables involved in the basic equations (0). The 
infinitesimal generator of such a group is a vector field X on Q 

whose components £ fc and r] a are supposed to be functions of class C°° on Q. The infinites- 
imal invariance criterion of the equations F = (F^F 2 ) under the group G is given by the 
determining equations 

(^ F Lo = ' ( 3 ) 

Y = t) {J) {W a )d uW (4) 
where means total differentiation 

D k = d k + u^du^ (5) 

and W a = r] a — ^ l u1 are the characteristics of the canonical vector field Y which is equivalent 
to X 

Y=X- £ k b k . (6) 

The determining equations Q are realized as an over-determined system of linear ho- 
mogeneous partial differential equations with respect to the unknown functions £ fc and rj a . 
In the determining equations the up to second order derivatives are the independent 
variables and these equations must be fulfilled identically over them on the manifold F. 

A sufficient condition of an existence of conservation laws is an invariance of the elemen- 
tary action associated with the Euler equations under a symmetry group of these equations. 
The starting point is the so-called Noether identity 1 

Y + Di ■ C = W a E a + biN\ (7) 
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In Eq.flJJ) N l are the Noether operators given by the expressions 

& = ? + D i{s) (W<*)(-l) r D ir) d uta)(r) (8) 

and 

E a = (-If D {s) d u « s) (9) 

are the Euler-Lagrange operators. The variational symmetries on the Euler-Lagrange equa- 
tions may be found via the Noether identity 

Y(L) + A • (CL) = (10) 

with the Lagrangian density A 



L = L(9, tp t , Ok, fk) = S (cos 9 - 1) (p t - ft 



1 n / \ 2 1 / -> \ 2 

sin 2 #(v^) +^(V6>) 



:n) 



.2 V T ) 2 

The dot in Eq.flZJ means the differentiation rule A ■ C {L) = LD t (f ) + fA (£)■ 
Substituting (JTTJl into (fTUj) one obtain 

Y(L) = W^L + A(^ 2 )^L + D k (W 1 )dg k L + D k (W 2 )d Vk L, (12) 

where 

d e L = —Ssm9(p t - ^ sin 26 (w) 2 i = $ ( c °s# - 1) , 

d dk L = -P0 k , d^L = -p sin 2 9ip k . (13) 
After a little manipulation one get the following equation for the unknowns £ fe and r] a 

r, 1 (-S sin 9ip t - | sin 25 ( w) "J 

,2\ ,„ A /'/■A ,„ A f<?fc 



+5 (cos6» - 1) A [V ) ~ <PtD t U ) ~ VkD t U 



-P9 k (A (v 1 ) - t D k (e) ~ OiDk 
-Psin 2 6<p k (A (V) - VtD k (e) ~ ¥>iD k (?)) 

+ (s (cos 9 - 1) cp t - I sin 2 9 (yvf - | (yfl) 2 ) (a (?) + A (?)) = 0. (14) 

Hereinafter, we keep the upper index k just for the space coordinates, the corresponding 
time component will be written explicitly. The coefficients at the derivatives of the functions 
9 and ip constitute an over-determined system of linear homogeneous partial differential 
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TABLE I: Equations for the uknowns £ fc and rf 1 . 



Independent 

V < I -L ±CAf KJ ±V_^ 


Equation 


o 


n? = 
'/t u 




f fn« /9 1^1 Tp 1 -\- rlinF 1 — rjl cnn $ 

^tUo U ) \ l(f ^ / — ' Dili 1/ 


CO I. 


-S 1 f cos 6* - 1) rf - 8 sin 2 6>77 2 = => rf = n? = 




'Ik u 




rf = r? 1 cot # + i ( rf + (iivfl + r? 2 


Muiflik ^ I) 


sin 2 6 ( rf + rf) = => rf + rf = 




e k = r,i + h(e t +divt) 


k Oi(k + I) 






^ = o 


(Vcp) (fk 


ej = 



equations with respect to the unknowns £ l and 77°. For example, extracting the terms at the 
independent variable 9 t and set equal them to zero 

-5 (ve - <Pk$) (cos9 - 1) + P0 k (e k + 6 k $ + ip k ($ 

+ (5 (-1 sin 2 6 (Vcp) 2 - \ (w)*) $ = (15) 
one obtain the conditions 

Ve =$ = £ = $, = $ = (16) 
that immediately imply £' = The results of the full analysis of the Eq. are 

summarized in Table 1. Finally, we get £' = £*(£), £ fe = 4^ fc (^)? 7 ? 1 = r ? 1 (^? 0)-, = ^(v 9 )- The 
uknowns £ fe satisfy to Killing's equation £j[ + £j = o - ^. As is well known, it has a solution 

for an arbitrary a just for n > 3. However, for the system considered there is a solution in 
2D case. Indeed, from the Table 1 one have 

<T = 2t* = 2nl + £t + divg, 

a = 2^ = 2T) 1 cot 9 + C + div£ + 2rf , 
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TABLE II: Conservation laws. 



Conservation 
law 



D, [ § sin 2 6 (Vy>) + f (w) j = A (W* + /?<^ fc sin 2 i 
A (S (cos 9- l)<p k ) = A (Tjh) 



A 



r x 



5 (cos 9 - 1) V^] J = A (e im "x m T nfc 
A (S (cos - 1)) = A (/?¥>fc sin 2 6) 



(l = z) 



(no summation over k) that imply 

e* = £*(*), m = s+Br, b t 



-B 




T] 1 = 0, if 



V. 



The solution found results in the following expressions for the characteristics W l = — £*# t 



a\j6 -b fx\/6 and W 2 = v - £ip t - cfcjip - b 



r x \jip 



conservation laws through Noether's theorem procedure 

AiV* (L) = 



. One can derive corresponding 



(17) 



under the condition E a L = 0. The Noether operators modify the Lagrangian density into 
a conserved quantity C l = N* (L) with a zero total divergence. Using the found variational 
symmetries we first construct the operators N l from formulae (JHJ) and then calculate from 
(JT7jl the corresponding conservation laws. The results are brought together in the Table 2. 
The physical interpretation of the obtained relations is obvious. 

1) The symmetry group of translations in time t — > t + a yields to the energy conservation 
law with the energy density W = sin 2 9 fVy>) +^f~ (^^) an d the energy density current 

f = -(js 2 e t ve + js 2 p t v^sm 2 e). 

2) This law connects the momentum P = hS (1 — cos#) Vy? and the canonical energy- 
momentum tensor 

T kl = hL5 kl + hp9 k 9i + h(3ip m sin 2 9 (18) 
and reflects the variational symmetry under space translations x\ — > x\ + a\. 



3) The conservation of angular momenta L 
tensor under the space rotation around z-axis. 



r x P 



through the angular momenta 



7 



4) The conservation law corresponding to the transformations tp ^ tp + \ connects the 
scalar variable N = S (1 — cos 9) with the current j = pVpsin 2 9. 

It is useful to apply the obtained laws to some types of solutions of Eqs. ((H). We restrict 
consideration by finite-amplitude spin-waves, dynamical soliton and skyrmion. 

1) The finite-amplitude spin- waves 9 = 9q = const and tp = kf — ufi with the dispersion 
huj = JSk 2 cos 9 . The last conservation law may be interpreted as a conservation of the 
magnon density N = S (1 — cos9 ). The corresponding magnon density current is j = 
^p-fcsin 2 #o- The space components of the energy-momentum tensor are 



kl 



S (1 — cos 9 ) fiui 
+ ifsin 2 # {kl-kl 

JS 2 sin 2 9ok x k v 



JS 2 sin 2 9 k x k y 

S (1 — cos 9q) Uuj 
-^sin 2 fl (kl ~ k 2 



The momentum P = S (1 — cos#o) hk is a times of the magnon density and the elemen- 
tary momentum hk. The energy density W = \ JS 2 k 2 sin 2 9 « JVftu; (0„ < 1) has the 
corresponding energy density current T = JS 2 ku sin 2 9 . We note a non- additive character 
of energy density at finite values 9 . 

2) For a dynamical soliton with axial symmetry one have to use the following parametriza- 
tion 9 = 9(r) and tp = The magnon density N = S (1 — cos#(r)) has the zero magnon 
density current j = 0. The energy-momentum tensor is 



ki. 



—S (1 — cos 9 (r)) ftu; 



+ ^f- cos 



(2urf) 



V 



JS 2 



M (f 



2 



-S 1 (1 — cos # (r)) 



-^cos(2^) (f)^ 

The momentum P is zero, the energy density W = and the energy density current 



T = 0. 



3) For the skyrmion the following parametrization 9 = 2 tan 1 (Jp) and tp = tan 1 (f) 7 
is used (for simplicity we consider an unit topological charge). Then, iV = 2S r2 ^ R2 , j r = 0, 



JS 2 



4R 2 r 



h (R 2 +r 2 ) 



The momentum in polar coordinates of xy-plane P r = 0, P v 



27i5 R 2 



however, T k i = 0. The energy density W = AJS 2 j. - 2 , 

-\-T ) 

current T — 0. 



TrR 2 



r r 2 +R 2 ' 

& and the energy density 
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III. DYNAMICAL INVARIANTS OF 2D TWO-SUBLATTICE MAGNET. 

In this section Lie transformation group methods will be applied to the partial differential 
equations describing the dynamics of 2D generic antiferromagnet (two-sublattice magnet). 
The established variational Lie symmetries for ferromagnet will be employed to derive anti- 
ferromagnet conservation laws revealing its important features. 

The action of 2D two-sublattice magnet is^ 

2 

L = S ( cos &i — 1) fit — P (—2 sin 9 1 sin 2 cos (<p\ — tp 2 ) — 2 cos 9\ cos 2 

+ cos 9\ cos 9 2 cos fa — (p 2 ) (V9\ V6* 2 ) — sin 9\ cos 9 2 sin (tpi — (p 2 ) (Vy?i V02^ 

+ cos 9i sin 9 2 sin ((p 1 — (p 2 ) (V9 1 V^J + sin X sin 9 2 cos ((pi — (p 2 ) ( Vy?i V^J 

+ sin 0i sin 2 (V0iV0 2 )) • (19) 

The time symmetry group t — » t + a under which the elementary action is invariant has 
the infinitesimal symmetry generator X = and Noether's operators iV* = 1 — @itw. — 

i=l lt 

2^2 2 

? ^**afc> Nk = ~ 5 ^«aSfc ~ ? ^*a£r- By introducm g the functions 
C* = AT* (L) = -/3 (-2 sin X sin 2 cos (^ - y? 2 ) - 2 cos 0i cos 2 

+ cos 0i cos 2 cos (</?i — ip 2 ) (V9i V0 2 ) — sin 0i cos 2 sin (ipi — y? 2 ) (V^i V0 2 
+ cos 0i sin 2 sin (ipi — ip 2 ) (V0i Vy? 2 ) + sin 0i sin 2 cos (fi — ip 2 ) \ VipiV(p 2 , 
+ sin 0! sin 2 (V0iV0 2 )) (20) 

and 

C7 fc = iV fc (L) = p (cos 0i cos 2 cos fa - ip 2 ) + sin 9 1 sin 2 ) (9 lt 9 2k + 2t 0i fc ) 
+(3 sin 0i sin 2 cos fa - <p 2 ) fa t <fi2k + f2t¥ik) 
+(3 cos 0i sin 2 sin fa - <p 2 ) (9 lt ip 2k + 9 lk ^2t) 

-P sin 0i cos 2 sin fa - ip 2 ) {9 2t (p lk + 9 2k (p lt ) (21) 
one get the energy conservation law A (C*) + D k [C k ) = 0. By the same manner one can 



obtain the generators and Noether's operators of the translation group x k — > x k + a k 

2 2 „ 2 2 

The momentum conservation law is 

£ A (5 (cos 0, -1) ^) = A (T w ) , (22) 
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with the tensor 



2 2 
T kl = L5 M - 9jhde jk L - ^ Vjkd^L. 

3=1 3=1 



The conservation of angular momenta is 



(23) 



]T A ([r x 5(cos^- - 1) V^] ) = D k (e znl x n T lk ) , 



3=1 



(24) 



X = e znl x n Sr 



N* = - e znl x n e\ 



dx l ' 

8 8 

— € X (ft 



d6f 



d d 

N k = znk n _ znl nga znl I n a 



The action is also invariant under a simultaneous rotation in both sublattices Lp\ = ip\ + x, 



2 2 - 

(p' 2 = y?2 + X- The generator and Noether operators are X = J2 — J2 X k = 

j=i (fj 



3 ~ x d<p jt ■ 



J2 jr~ and they imply the conservation law 

j — l a ¥jk 



^2 D t (S cos 9j) = (3D k (sin 9 1 sin 6 2 cos (<pi - <p 2 ) (<fik + <-P2k) 

3=1 

— sin (</?i — ip 2 ) (9 2k sin 9\ cos 9 2 — 9± k sin 9 2 cos #i)) . (25) 



IV. CONCLUSIONS. 



Lie transformation group methods have been applied to the partial differential equations 
describing the dynamics of ferro and antiferromagnets and the variational symmetries of 
the action functional on the Euler-Lagrange equations are established. Each such symmetry 
gives rise to a conservation law with a conserved current. The explicit expressions for the 
conserved currents are found. 



Appendix. 

Let's take the following notations: 

1) the matrix of Frechet derivatives fg — ^2 -^-D( s ), where the sum is taken over all 

S " U ( s) 

multi-indices s — (si, s 2 , . . . , s k ) with s k — 1, . . . , 3 (k > 0); 
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2) the adjoint operator of Frechet derivatives- 2 * 



df a 



One have to found a Lagrangian density L (x,u,u , . . that the initial set of differential 
equations are the Euler-Lagrange equations of the following action 



t-2 



From the ratio 



E a L = r+ f-f 



Xu 



x 



dX 



is seen that that sufficient condition of the solution of the problem is the self-adjoint condition 
of the /-operator 



df a 



-l) s+r C°D 



OU (s)(r) 



or in the matrix notations f@ = (/^ = • The straight calculation yields to 



./ 



-Sip t cos 6 - (3 (Vtp) 2 cos 26 + (3 A -S sin 6d t - (3 sin 26 (V^V) 



S6 t cos + 2/? cos 26 (Vy>V0) + /? sin 26Atp 
+S sin 0d t + /3sin20 (V^V) 



/3 sin 20 (v^v) + (3 sin 2 A 



and 



-Sfo cos - 13 (y<pf cos 20 + (3 A 
\ -S sin 6d t - /3 sin20(v^V) 



S0 4 cos + 2/3 cos 20 ( Vy?V0 + /3 sin 20A</? 



+5 sin 6d t + (3sin26 (V</?V 
/3 sin 26 (V^V) + /3sin 2 0A 



/ 



Obviously, the claimed condition is fulfilled. The Lagrangian density can be constructed via 
the homotopy formula 2 - 

l 

L[ u ] = ]>> a (z) / f a (x, Xu, Xu',.. .) dX. 
o 

Using the explicit form of differential equations it results in 

i i 
L = —S (dipt - <p9 t ) J dXX sin X6 + [36 Ad J dX X 
o o 

+13 (<p (V^W) - - (V^) 2 ") J dX A 2 sin 2X6 + (3<pAip J dXX sin 2 X6, 
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or after some manipulations in 



L = S(cos0-l)</> t -|sin 2 0(v^) 2 -|(W) 



+SD t 



sin 9 \ 



where 



C = ^[6 (6 - sin 20) + sin 2 9 



The last two terms present total derivatives and may be dropped. 
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